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Rephrasing the backbone of two-dimensional percolation as a monochromatic path crossing 
problem, we investigate the latter by a transfer matrix approach. Conformal invariance 
links the backbone dimension to the highest eigenvalue of the transfer matrix T, and 
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1. Introduction 

The critical behavior of percolation has attracted considerable interest in the mathe- 
matical physics literature over the last decades. Whereas most practical applications (such 
as studying the efficiency of oil extraction from a porous soil, or the fractal geometry of a 
strike of lightning) take place in three spatial dimensions, analytical progress has largely 
been confined to two dimensions Although of geometric origin, percolation fits in 
the framework of critical phenomena, and in particular the concept of universality should 
apply. One therefore expects the specific choice of a discrete model (bond or site percola- 
tion) and of the lattice structure (e.g. square or triangular) to be of no relevance to the 
determination of the critical exponents. 

A large part of the progress made is due to the identification with the g ^ 1 limit of 
the g-state Potts model [|l|] . A very fruitful idea has been to treat the latter in terms of its 
random cluster formulation 0, and further in terms of the loops surrounding the clusters 
0]. Applying Coulomb gas (and related) methods to the loop model led to a range of exact 
results around 1980 |^. In particular, the correlation length exponent = | and the 
magnetic exponent Xh = ^ MM (^^^ codimension of which is the fractal dimension of the 
percolating cluster, D = 2 — x/^ = ||) were computed. 

The next major advance followed from the advent of conformal field theory 0], which 
provides an appealing correspondence between the g-state Potts model (for particular 
values of q) and the so-called minimal models. For instance, the exponents Xk = -^{k"^ — 1) 
with k > 2 [PIPIH] describing the asymptotic decay of the probability Pk{r) ~ r"^^*" of 
having k loop segments connecting two narrow regions over a distance r ^ 1 were 
found to fit in the Kac table of conformal dimensions ITlIl. Another remarkable result is 



the celebrated Cardy formula |T^ expressing certain path-crossing probabilities in terms 
of hypergeometric functions. 

More recently, percolation has attracted the interest of probabilists. In a ground- 
breaking publication, Smirnov has proved that the scaling limit of site percolation on the 
triangular lattice exists and is described by the stochastic Loewner evolution with param- 



eter K = 6 Consequently, most of the results referred to in the above have now been 
rigorously proved. 

Nevertheless, a certain class of exponents have continued to resist the physicists' at- 
tempts over the years. These are most conveniently defined by considering bond percolation 
inside a large square, of which we imagine two opposing sides to be connected to super- 
conducting plates (see Fig. 1). Each percolating bond is stipulated to possess a fixed and 
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finite conductivity, and an electric voltage is applied across the plates. At the percolation 
threshold p = pc, the part of the network that supports a non-zero current is known as the 
backbone, and its fractal dimension determines a critical exponent X2 = 2 — L>b- Near 
Pc, the conductivity of the network scales as {p — PcY, defining the conductivity exponent 
t. The latter can be connected to the fractal dimension of random walks constrained to 



the percolating cluster, or to its backbone, via the Einstein relation [|I4 . 

A number of conjectures for X2 have been falsified as numerical simulations have 
become increasingly accurate. The benchmark thus far is the Monte Carlo method of 



Grassberger [[151 in which the conducting part of the cluster is identified using a clever 
recursive algorithm. Large-scale simulations yield the value X2 = 0.3568 ±0.0008 The 
exponent X2 is actually a member of a family of so-called monochromatic path-crossing 



exponents Xk [rD||, with the magnetic exponent fitting in as Xh = xi. The higher exponents 



Xk, k > 3 are all unknown. 

In the present publication we provide a numerical estimate of X2 using an algorithm 
which is entirely different from that of Grassberger. Using the reformulation of X2 as a 
path-crossing problem, we relate it to the largest eigenvalue of a linear operator (actually 
a transfer matrix) that builds all possible percolation clusters supporting at least k = 2 
mutually non-intersecting paths. We work in the geometry of semi-infinite strips of width 
L, with L < 9. 

Our approach is interesting in several respects. First, the reformulation as an eigen- 
value problem makes direct contact with the predictions of conformal field theory ||17|| . 



That Grassberger 's recursive algorithm defines a conformally invariant observable is not a 
priori obvious, but the fact that the transformation to a path-crossing problem involves a 
conformal transformation and that we here obtain a consistent value of X2 shows that this 
is indeed the case. One would then further expect X2 to be the conformal dimension of a 
primary operator O2 in some (presently unknown) conformal field theory of percolation. In 
particular, the conformal tower of O2 should possess descendents whose conformal dimen- 
sions are integer-spaced with respect to X2- We have checked this prediction by examining 
the scaling of the first few eigenvalues of our transfer matrix with system size. We shall 
present evidence of a level two descendent with conformal dimension 2.35 ± 0.1, whereas 
there does not appear to be a descendent at level one. 

Second, the generalization of our method to the case of more [k > 3) paths, or to 
the Potts model with q 1 states, are immediate. Results for these cases will appear in 
a separate publication [^. Third, from a technical point of view we have had to tackle 



the major obstacle of writing a transfer matrix in which some degrees of freedom (the 
percolation clusters) must be summed over, whereas others (the paths) act as constraints 
on the former but must not themselves be summed over. Fourth, we have devised an 
algorithm which is naturally parallelizable. 

Like Grassberger we find that the data for X2 are hampered by strong (presumably 
non-analytic) corrections to scaling. As a consequence our final result 

;r2 = 0.3569 ± 0.0006 (1.1) 



confirms that of [|T6[ , but unfortunately does not improve its precision. On the other hand, 
we have devised some variants of our algorithm in which the constraint of mutual avoid- 
ance of the two paths is relaxed, so that they are allowed to touch in some configurations 
at vertices but not to share an edge. Physically this means that we measure the fractal 
dimension of the backbone with some "blobs" that are tangent to it included. The sur- 
prising result is that this relaxation of the original definition does not alter the value of 

The paper is laid out as follows. In Sec. 2 we review the reasoning leading from the 
original formulation of the backbone dimension to that of a path-crossing problem, and 
we restate the latter in a strip geometry. The construction of the corresponding transfer 
matrix, and of its associated state space, is described in Sec. 3. In Sec. 4 we transcribe this 
as an algorithm and discuss its implementation. The data is analyzed and extrapolated 
to the L —>■ oo limit in Sec. 5. The appendix displays some transfer matrices produced for 
small system size. 

Note added: When this work was being completed we became aware of the preprint 
of Lawler, Schramm and Werner in which iii = ^ is established on a rigorous basis, 
following Smirnov [|TB[. The authors also relate X2 to a second-order partial differential 
equation with specific boundary conditions, but fail to provide an explicit solution of the 
latter. We thank John Cardy for bringing this to our attention. 



2. Path-crossing probabilities 

Let us return to the formulation of the backbone problem given in the introduction, 
namely in the so-called busbar geometry (see Fig. 1). The condition that a given point 
(site or bond, as the case may be) on the spanning cluster belongs to the backbone is 
that it can be connected to either of the superconducting plates by means of two mutually 
non-intersecting paths. 

^ Strictly speaking, this condition includes also points which are being held exactly at zero 
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Fig. 1: Busbar geometry, here shown for the case of bond percolation on the 
square lattice. The backbone is indicated by fat edges. 

This choice of geometry is somewhat unnatural, as it does not fully display the ro- 
tational symmetry of the continuum limit. It is more convenient to work in an annular 
geometry limited by two concentric circles of radii r ^ 1 and R ^ 1. Interpreting the 
inner circle as the point which is a potential element of the backbone, and the outer cir- 
cle as the point at infinity, we see that a given percolating configuration in the annulus 
contributes to the backbone if and only if the two circles are connected by two mutually 
non-intersecting paths on the percolating cluster (s); see Fig. 2. 

More generally, one may define higher exponents Xk hy studying, at the percolation 
threshold, the probability Pk{r, R) ~ (^)^'" that the annulus is traversed by k mutually 
non-intersecting paths. Clearly, the configurations in which these paths belong to different 
clusters are asymptotically sub dominant, and so we might as well assume that they belong 
to the same cluster. 

The situation may be further refined by considering path-crossing events in which 
a given number of traversing paths belong to the clusters, and the remaining number belong 
to the dual clusters .11 More precisely, for each k there are 2^ types of path configurations, 

current in a Wheatstone's bridge- like arrangement. Since in the continuum limit the percolation 
cluster is almost surely "asymmetric", such points are extremely rare. See also the discussion 
below on the possibility of contact points for paths. 

^ For site percolation, the dual clusters consist of the non-conducting (uncolored, white) sites. 
For bond percolation on the square lattice, it is most natural to think of the dual clusters in terms 
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Fig. 2: Annular geometry endowed with critical percolation (here shown in 
the continuum limit). A possible choice of two disjoint percolating paths is 
shown as dashed lines. 



each specified by a set of color variables (ri, r2, . . . , Tfc) with = +1 (resp. = — 1) 
meaning that path number i belongs to the clusters (resp. to the dual clusters). Within 
the context of the g-state Potts model (with q ^ 1), it is not obvious whether different 
choices of the color variables will lead to the same critical exponents, except of course for the 
obvious symmetries obtained by rotating the sequence (ri, r2, . . . , r/j), reversing its order, 
or dualizing it. But in the percolation case {q = 1) the bonds (or sites) are uncorrelated, 
and various parts of the system may be dualized independently. Using this approach, it 



has been proven in the case of site percolation on the triangular lattice that all the 
polychromatic sequences (in which both a = +1 and a tj = —1 are represented; k > 2) 
share the same critical exponents. In particular, any polychromatic color configuration 
may be transformed into the alternating one, = (—1)*. 

We expect this result to be independent of a particular lattice realisation, and thus 
to apply also to bond percolation. In this case, the identification of the critical exponent 
with that of k traversing loop segments on the surrounding lattice, referred to as Xk in 
the introduction, becomes evident (at least for k even). A rigorous proof that the formula 
Xk = -^{k"^ — 1) applies to the polychromatic path crossing problem for site percolation on 
the triangular lattice was spelled out in [jl^ . 



of the standard duality transformation in the random cluster model according to which any 
conducting edge is intersected by a non-conducting dual edge, and vice versa. 



5 



For monochromatic sequences (all = +1) the argument given in [10] fails, and the 
corresponding exponents Xk are expected to be different from the Xk- Indeed, from entropic 
considerations it should be clear that Xk < Xk < X2k- 

Several of the Xk have nice physical interpretations. Thus, X2, x^ and 0:4 are re- 
spectively the codimensions of the cluster perimeter (hull) 0, of the external (accessible) 



perimeter , and of the set of pivotal (singly connecting) bonds . The latter also yields 



the correlation length exponent 0, via the scaling relation u = 1/(2 — X4). 

In the absense of an exact solution, one might imagine evaluating the exponents Xk 
numerically by measuring the decay of the path crossing probabilities on an annulus, as 
outlined above. A more feasible alternative is to compute certain restricted free energies on 
semi-infinite cylinders by using a transfer matrix, as we shall describe in the next section. 
These free energies can be related to the critical exponents as follows. 

Since the scaling limit of critical percolation is conformally invariant P,13|, one is 
allowed to transform the annular geometry of Fig. 2 into a cylindrical one by means of the 
conformal mapping w = u + iv = log{z). The transformed complex coordinate w may 
be thought of as imbedded in the strip —00 < u < 00, < f < L with periodic boundary 
conditions in the f-direction. All this means is that Fig. 2 must be viewed in perspective, 
interpreting the inner and outer circles as the extremities of the cylinder. 

We are going to make use of the following result: let fo{L) be the free energy per unit 
area for the unrestricted percolation problem, and fk{L) (resp. fk{L)) the corresponding 
quantity for the constrained problem where only those configurations are included in the 
partition sum in which (at least) k monochromatic (resp. polychromatic) paths span the 
length of a semi-infinite cylinder of width L. Then as L — > 00, the discrete lattice model, 
at criticality, should have a continuum limit described by conformal field theory, so that 

m 



A(L)-/o(L) = ^ + o(L-2), (2.1a) 
/fc(L)-/o(L) = ^ + o(L-2). (2.16) 



We shall obtain estimates for the Xk hj extrapolating data for sufficiently large strips 
to the limit L — 00. 
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3. Transfer matrix algorithm 



It has been known for a fong time how to numericaUy compute the fk{L), by writing 
the transfer matrix for the loop model in the basis of planar (Catalan-like) connectivities 
(see 1^ for a closely related computation). The same is true for fi{L) by using the trick 



of adding a ghost site |j2T|], or alternatively (via a duality argument) by forbidding the 
clusters to wrap around the cylinder . 

The computation of f2{L), the principle of which we now describe, is considerably 
more complicated. The main complication stems from the fact that to compute the cor- 
responding partition sum we must exclude those configurations of the percolation clusters 
that do not support (at least) two spanning paths, and count each of those that do with 
unit weight (and not with a weight equal to the number of ways two such paths can be 
realized for the given cluster configuration). Roughly speaking, the degrees of freedom are 
the clusters and the paths, and we must trace over the former but not the latter. 

For the sake of definiteness we consider in this section critical bond percolation on a 
square lattice, though the principle of the transfer matrix can be applied to any lattice 
with any probability of occupation p, and to bond as well as site percolation. Since in our 
case Pc = ^ Q] ' convenient to simply assign a weight of one to every configuration of 
percolating/non-percolating bonds. For now we consider the simplest orientation of the 
lattice, which corresponds to L sites in the transverse direction with periodic boundary 
conditions (Fig. 3). With all these conventions, fo{L) — —2 log 2 in Eq. (2J.). We shall 
discuss later another possible orientation of the lattice. 
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Fig. 3: The square lattice with periodic boundary conditions along one of 
its orientations. The dotted lines are time slices. 
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Fig. 4: The 10 possible path configurations for L = 4. 

We keep track of the paths by defining path configurations in analogy with those used 
in the transfer matrix calculations of the self-avoiding walk among the L sites in a 



row (at time t = to), two sites are connected to the point at infinity (time t = — oo) by 
means of paths. Furthermore, in order to allow subsequent backtracking of either path 
(at a later instant t > to), the remaining sites may be connected in pairs by means of 
backward arches. The possible configurations for L = 4 are listed on Fig. 4. 

To overcome the difficulty of not summing over the possible path configurations, we 
define the basis states on which the transfer matrix acts as lists of path configurations. 
Elements of the list give all possible realizations of the positions of the paths (and of the 
arches) which are compatible with the "past" of the state. 

Formally, if V is the set of path configurations, then basis states are indexed by non- 
empty subsets of V (one must exclude the empty subset since it corresponds to states for 
which there is no possibility of two disjoint paths reaching time t) . Note that the dimension 
of the total space is — 1, which grows extremely rapidly with L. We shall return to 
this point when we discuss practical implementation. 

By definition, the matrix element T^B between basis states indexed hy A G V and 
B <Z V equals the number of configurations of the bonds between time t and time t + 1, 
such that the state A at time i -fl is obtained from the state B at time t. Given the initial 
state B and the configuration of the bonds a; G O (that is whether they are percolating or 
not, O being the set of all possibilities), the procedure to determine the final state A is as 
follows: 

★ For each possible initial path configuration h E B, consider all possible continuations of 
the existing lines at time t (the two original paths and the arches) that are compatible 
with the configuration of the bonds uj. Note that each line must be either continued 
to a site at time t or be connected to another line (in which case it will reemerge 
at the other end of the arch; the lines coming from infinity or from the same arch 
cannot be connected to each other). Furthermore, for each pair of adjacent empty 
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sites one must consider the possibility of creating a new arch. Let 0(6, uj) dV he the 
hst of path configurations at t+1 thus produced. 

The full state A is reconstructed by simply putting together all the possibilities (of the 
form (j){h,ijj), b e B) obtained for each initial path configuration. If one finds A= 
this means that no continuation is possible, and the state is excluded. 
We give an example of such a computation on Fig. 5. 



{ 




{ 



1^ 



T ) 



Fig. 5: Evolution of two path configurations with the same percolation 
configuration. Solid (resp. dashed) lines represent percolating (resp. non- 
percolating) bonds, whereas thick lines represent the possible paths. 

In other words, we have the formal identity: 

Loeii beB 

which shows quite explicitly that one sums over bond configurations but not over path 
configurations. 

Finally, the free energy per unit area is given by 



/2(L) = - lim -log(^|T(L)*|B> 

t — >00 ±jI, 



(3.1) 



where the states A, B specify the boundary conditions and are essentially arbitrary (the 
state A should belong to the image of T, see next section), and T(L) is the transfer matrix 
for strip width L. 
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As t — > oo, the matrix element T(L)* is dominated by the largest eigenvalue 
A(L) of T(L), and combining Eqs. ( |2.1a| ) and ( |3.1| ), we find: 

1 A(L) = l-^ + o(L-i) (3.2) 
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4. Algorithmic details 

In order to appreciate how effective the transfer matrix approch is, it is important 
to understand the structure of the matrix constructed in the previous section. It is an 
integer-valued matrix of extremely large size, but many of its entries are zero. In fact, 
starting from any basis state a very limited number of states are generated. These 
are the only states that matter for the determination of the largest eigenvalue(s) and we 
can thus restrict ourselves to a submatrix of much smaller size. 

We now describe schematically the procedure we used. The main steps of the algorithm 
are as follows: 

(i) Start with an arbitrary basis state (ideally, one that we know is generated by iteration 
of the transfer matrix). Put it onto a "stack" of states to process. 

(ii) Pick a state B from the stack and "process" it, i.e. generate the non-zero entries Tab, 
and store them. This encodes one column of the transfer matrix. 

(iii) Consider every new basis state A that has been generated at step (ii); check if it has 
already been processed; if not, add it to the stack. If the stack is non-empty, go back 
to step (ii). 

(iv) Finally, once the stack is empty, the largest eigenvalue is computed by simple iteration 
of the matrix that has been generated. 

The transfer matrix is such that the submatrix thus generated has no zero rows or 
columns. We call this submatrix the reduced transfer matrix. 

An important remark for practical applications is that this procedure is highly paral- 
lelizable: several CPUs can perform step (ii) simultaneously and independently, only the 
stack must be shared. In practice, it is necessary to have a server that communicates with 
the various clients involved in the computation; it ensures that their stacks are synchro- 
nized, and dispatches the tasks. At the end of each calculation (step (ii)), a client sends 
the server the new states created and receives the states created by other clients in the 
meantime. The time spent updating the stack being very small compared to the calcula- 
tion time, the parallelization is near 100% efficient (at least up to 20 clients which is the 
maximum we tested). 
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Let us now discuss in more detail this procedure. 

First, we must define how to encode path configurations. A study of Fig. 4 shows 
that if exactly k = 2 paths are connected to t = — oo, then they can be considered as 
an extra arch. This trick reduces the number of configurations and slightly simplifies the 
implementation (but cannot be extended to /c 7^ 2). We can then move the point at infinity 
and redraw the configurations as standard arch configuration^, see Fig. 6. 



Fig. 6: The configurations of Fig. 4 redrawn as arches. 

An arch configuration is then encoded in a standard way as a sequence of closing / 
empty / opening steps, that is G {—1,0, +1}, I < i < L, such that the height function 
hi = X]f=i satisfies hi > for all i and /il = 0. States are now defined as sorted lists of 
path configurations. 

Next, we discuss how to perform step (ii) in practice. One possibility would be to 
apply directly the principle of section 3, that is to consider all possible bond configurations 
between 2 successive time slices and for each, to produce the resulting state. However, 
since there are 2^^ such configurations, the time required to do so grows exponentially, 
which is not satisfactory. Besides, the determination of all possible continuations of the 
paths to time t + 1 is a rather complex task. Instead, we shall use a factorization of the 
transfer matrix as a product of L sparse matrices T^, 1 < i < L which describe the addition 
of a single site. The details of the factorization depend on the exact situation envisioned. 
We present here three cases. 

4.1. The square lattice with standard orientation 

The example used so far is that of the square lattice with its usual orientation. In this 
case the factorization can be pushed further by writing that T(L) = Hi . . . Vi . . . 
where (resp. H^) corresponds to the addition of a single vertical (resp. horizontal) bond 
(Fig. 7). 

^ The number of L-point arch configurations equals m_L — 1, where rriL are the Motzkin numbers 
(the empty configuration is excluded). The generating function M{x) = XliT-o'^^'^^ ~ 
(1 — X — \/l — 2x — 3x'^)/2x^, has a singularity in x = Xc = |, showing that the number of path 
configurations is ~ 3^ asymptotically. 
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2 H=Hi H2 ...H 
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J V=Vi V2 ...Vl 
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L 



Fig. 7 : Factorization of the transfer matrix. 

The action of Vj is very simple: Vj = + V" where (resp. ) describes the 
evolution when the vertical bond number i is percolating (resp. non percolating). is 
simply the identity, whereas Vf acts on path configurations as follows: either a path/arch 
is at site z, in which case it gives (the path cannot cross the non-percolating bond), or 
there is not and it is the identity. The action on a state made of several path configurations 
can be deduced from these basic rules, as explained in section 3. 

The action of is slightly more complicated: Hj = + H", similarly as above. 
The definitions of and H." must take into account all the possibilites of continuations, 
recombinations and creations of paths along the horizontal bonds. This requires working, 
as intermediate states, with path configurations of length L + 2 instead of L, since one 
must temporarily distinguish the paths directed horizontally and vertically at the first and 
last vertices being currently processed. We leave the details as an exercise to the interested 
reader. 

A final ingredient is that one can use the dihedral symmetry of the transfer matrix: 
since the latter commutes with cyclic permutations of the sites and with refiections, one can 
select a representative in each orbit of the dihedral group among the basis states. Note that 
the action is an overall action on all configurations that constitute the state simultaneously. 

The states generated by the procedure above can then be replaced with the representative 
state of their orbit, producing a smaller transfer matrix but with identical eigenvalues. 
This further reduces the size of the transfer matrix, by a factor of (roughly) L. 

4-2. The square lattice with standard orientation 2: the square/octogon deformation 

It is interesting to study variants of the algorithm above. One natural question is: if 
one allows the paths to touch each other at vertices, how is the asymptotic behavior of the 
free energy modified and in particular is £2 left unchanged? Another possible formulation 
of this question is to consider a deformation of the lattice in which each vertex is replaced 
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with a small square, resulting in a square/octogon lattice (Fig. 8). The bonds of the small 
square are always percolating and allow paths that would have touched at a vertex to avoid 
each other E 




Fig. 8: Deformation of the square lattice. 

Physical insight suggests that such modifications should not affect the values of X2- 
The reason is that, just like the Wheatstone bridges configurations mentioned in the in- 
troduction, the fact that current fiows through loops which are connected to the backbone 
by just one point is rather unstable since any microscopic defect that breaks the symme- 
try between the two orientations of the loops (deforming the lattice is precisely a way of 
introducing such a defect) will produces a non-zero current. If X2 is to be universal it 
should not depend on such microscopic details. It is this insight that we would like to test. 
There is another, more practical reason one would want to study such modifications of the 
algorithm, which will be apparent in section 5. 

It is very simple to modify the transfer matrix of section 4.3 to allow such path 
evolutions. Vi is unchanged, whereas Hi now allows two paths to reach the same vertex 
and emerge from it as if they had not touched each other. 



4-3. The square lattice with light-cone orientation: the hexagon deformation 

Finally, we rotate the lattice by 45 degrees, the motivation being that we expect 



better convergence properties, as observed empirically in similar computations |22,25 



Unfortunately, there is no efficient way to encode the corresponding configurations, and 
we are therefore led to a modification of the lattice which is similar to what was done in 
section 4.2: this time the result is a hexagon lattice in which vertical bonds are always 
percolating (Fig. 9). This is equivalent to allowing "horizontal tangencies" on the original 



^ Note that a path crossing a vertex of the original lattice can correspond to two different paths 
on the deformed lattice, but since we do not sum over path realizations this is of no consequence. 
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Fig. 9: Another deformation of the square lattice. 

square lattice, that is allowing two paths to touch at one vertex in the configuration where 
the two upper edges belong to the same path; however, "vertical tangencies" are still 
excluded. 

In this case, encoding the states becomes completely identical to what was done pre- 
viously. There is a decomposition T = Ti . . . where adds an extra vertex i at time 
t + 1 (and two bonds). Since the new sites at t + 1 are now shifted with respect to the 
sites at t, the action of the transfer matrix includes a conventional rotation of a half-bond's 
length (or tt/L). 

Relations (|3.1| )-(|3?^) must also be modified to take into account the 45 degrees rota- 
tion; the latter introduces an extra factor of 2 in the unit of area, so that fo{L) = —4 log 2 
and: 

^A(L) = l-^+o(L-i) (4.1) 

This factor of 2 alone increases the accuracy of the measurement of X2 compared to the 
other two cases, since the corrections are expected to be smaller. 



5. Numerical results 

We show on table 1 the size of the reduced transfer matrix for 4 < L < 9, in the three 
cases presented above (sections 4.1, 4.2, 4.3). While the full matrix is very roughly of size 
2'^ , the size of the reduced matrix seems to grow as L\, which is still large but not as 
intractable. It is interesting to note that S2 < si, that is the modification of the lattice to 
allow configurations where paths touch at a point decreases the number of states. 



L 


4 


5 


6 


7 


8 


9 


Sl 


15 


72 


515 


4219 


41728 


? 


S2 


12 


51 


291 


1893 


14923 


132799 


S3 


12 


51 


291 


1893 


14923 


132799 



Tab. 1: Size of the reduced transfer matrix. 
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We have no deep explanation for the remarkable equality of sizes of algorithms 2 and 
3, except the observed fact that the states generated are the same in the two cases. 

Next we present the data for the largest eigenvalue of the transfer matrix on table 2 
with a twelve digit accuracy. 



L 


4 


5 


6 


7 


8 


9 


Ai/22^ 

A2/22^ 

A^/22^ 


0.514287790945 
0.540388840500 
0.718747415570 
0.058692638251 


0.594678112301 
0.617254658842 
0.775012703547 
0.145046191784 


0.653760363032 
0.672285202673 
0.812529692986 
0.224345992159 


0.698459489246 
0.713573950794 
0.839330907375 
0.292806902950 


0.733243927216 
0.745682316102 
0.859432882632 
0.351338353673 


? 

0.771356857232 
0.875067710677 
0.40153182 



Tab. 2: Largest eigenvalue of the transfer matrix. The last row shows the 
second real eigenvalue for the third transfer matrix. 

In order to study the asymptotic behavior of these series of numbers, we use Eq. ( |37^ ) 
for cases 1 and 2 (or (|4.1|) for case 3) to extract approximate values of X2- The results are 
on Fig. 10. We also presented quadratic fits of these data. 




Fig. 10: Values of X2 obtained from the eigenvalues of the transfer matrices 
(Tab. 2). The order of colors is: 1 - red, 2 - green, 3 - blue. 

Several remarks are in order. First the two curves corresponding to the square lattice 
with its regular orientation (with or without contacts at points) seem to converge nicely 
within the range allowed by the fits. This means that the value of X2 is not affected by this 
modification. However, it is clear that the next corrections to Ai and A2 are quite different. 
Secondly, it is again manifest on figure 10 that the third set of data, corresponding to the 
45 degrees rotated square lattice, reaches its limit much faster than the other two. Whereas 
various fits will give a limiting value for the first two anywhere between 0.355 and 0.36, 
the range is limited to 0.3563 to 0.3575 for the latter. Assuming all these limits to be the 
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same, we reach the estimate ( p. . 1|) mentioned in the introduction. Note that there is no 
simple way for us to evaluate error bars since the results are entirely dependent on the fits 
used, the latter being arbitrary without any knowledge about the subleading corrections. 

Finally, numerical estimates of the norms of higher eigenvalues of the transfer matrix 
spectra can be extracted by a standard iteration/orthogonalization procedure ||2^. Using 



this method, complex eigenvalues are characterized by an oscillatory behavior and can thus 
be discarded (we expect physical observables to be linked to real eigenvalues) . Specializing 
to case 3 (cf. section 4.3 above), we find the fourth eigenvalue (in norm) to be the second 
real one. Its finite-size scaling is well fitted by (|2.1|) , defining a critical index 



4 = 2.35 ±0.1 (5.1) 



This is consistent with the conformal dimension of a level two descendent of the backbone 
operator. 

Extracting the scaling dimensions for even higher eigenvalues becomes increasingly 
problematic, as the finite-size effects get considerably stronger. It should however be 
noticed that the third real eigenvalue is doubly degenerate for any width L > 4. This is 
supposed to have implications for the organization of the conformal tower of the backbone 
operator. 
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Appendix A. Structure of some small size transfer matrices 

As an illustration of the algorithm explained in this article, we provide here the sim- 
plest non-trivial transfer matrices obtained with the geometries of sections 4.2 and 4.3. 
They correspond to a strip length L = 4 and their size is s = 12. 
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Fig. 11: Basis states (up to overall dihedral transformations) for L = 4. 



On Fig. 11 is described the basis in which these matrices are expressed. 



The matrices themselves read: 
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18 6 6 12 8 
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6 26 30 32 14 11 
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